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Abstract 

Oh I This paper is devoted to hyperbolic systems of balance laws with 

■ non local source terms. The existence, uniqueness and Lipschitz de- 

fH | pendence proved here comprise previous results in the literature and 

can be applied to physical models, such as Euler system for a radiating 
gas and Rosenau regularization of the Chapman-Enskog expansion. 
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£2 ■ 1 Introduction 

O ■ 

This paper is devoted to systems of conservation laws with non local sources, 
^ . i.e. to equations of the form 

d t u + d x f{u) = G(u) (1.1) 

>'• 

• i-h . where / is the flow of a nonlinear hyperbolic system of conservation laws and 

G: L 1 \—* L 1 is a (possibly) non local operator. As examples, we consider 
below the case G(u) = g(u)+Q*u that enters a classical radiating gas model, 



C3 



see [21], as well as Rosenau regularization of Chapman-Enskog expansion of 
the Boltzmann equation, see [20l [21] . We establish that (jl.ip is well posed 
in L 1 , locally in time, for data having sufficiently small total variation. To 
this aim, we require on (jl.ip those assumptions that separately guarantee 
the well posedness of the convective part 

d t u + d x f(u) = (1.2) 

and of the source part 

d t u = G{u) . (1.3) 
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These two equations generate two semigroups of solutions, say S and S. 
To obtain our results we exploit the techniques in [2j [9] , essentially based 
on the fractional step algorithm, see [9l [TUl [I3j [23] . Its core idea is to get 
a solution of the original equation as a limit of approximations obtained 
suitably merging S and S. 

On the two semigroups we require the following two key conditions 
(see P (S2) and (C*)]): 

i) A Gronwall type estimate with respect to a suitable metric d(-, •): i.e. for 
a positive C 



d(S t u, S t v) < e ct d(u, v) 
d(E t u,E t v) < e ct d{u,v) 



ct „ nA for all t > 0. (1.4) 



ii) A commutativity relation 

d(E t S t u, StS t u) < K t 2 as t -> 0. (1.5) 

The former assumption is used to prove the uniformly continuous depen- 
dence of the approximations from the initial data, also called stability con- 
dition in the framework of Lie- Trotter formula, see [14} Corollary 5.8, Chap- 
ter 3]. The latter condition yields the convergence of the approximations 
and ensures the uniqueness of the limit. 

Assume that i) holds and d is a reasonable (in the sense of Proposi- 
tion [3]T8]) metric equivalent to the L 1 distance. Then, the invariance under 
the hyperbolic rescaling (t, x) — ► (Xt, Xx) of solutions to system of conser- 
vation laws, implies that C = 0, see Proposition 13.181 Hence, to apply the 
operator splitting techniques, we need a contractive metric for the conser- 
vation law (|1.2|) . This role is naturally played by the well known functional 
$ in [18} 119]. Note, however, that this functional is not a metric, for it may 
lack to satisfy the triangle inequality. The proof, then, consists in showing 
that the semigroup generated by the source part (II. 3|) satisfies (ll.4p with re- 
spect to Liu & Yang functional and commutes with the semigroup generated 
by the conservation law in the sense of (|1.5p . 



More precisely, let f2 be a open subset of M. n with £ fi. For all positive 
5, define 

Wj = {«e L^RjJtyTV^) < 5\ . 

As a general reference on conservation laws we refer to [H] . On the convective 
and on the source parts we assume throughout that 

(F) / E C 4 (S7;IR n ) is strictly hyperbolic and each characteristic field is 
either genuinely nonlinear or linearly degenerate. 
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(G) For a positive S , G:U$ i— > L 1 (IR, M n ) is such that for suitable positive 

\/u,w£Us a \\G(u) — G(w)\\- L1 < L\ ■ \\u — to|| L i 
VttGW Jo TV(G(n)) < L 2 -TV(u) + L 3 . 

Note that (F), respectively (G), ensures the local in time well posedness 
of (II, 2p . respectively (|1.3jl . A class of functions satisfying (G) is provided 
by the following proposition. 

Proposition 1.1 Let g,h:VL^ W 1 be locally Lipschitz and Q £ L 1 (M; M nxn ). 
Then, the operator G(u) = g(u) + Q * h(u) satisfies (G) with L3 = 0. 

The proof is deferred to Section [3j We are now ready to state the main 
result of this work. 

Theorem 1.2 Let f satisfy (F) and G satisfy (G). Then, there exist pos- 
itive T, 8, C, closed domains T>t and processes 

F t :V T - t ^V T Vte[0,T] 

with the properties: 

(1 ) for t,s £ [0, T] with t < s, U~ s QV t CP, C U 5o ; 

(2) for u in V T , F u = u; for t,s £ [0, T] with t + s £ [0, T], F s V t C V t+S 
and for u G V T ^ t - s , F t F s u = F t+S u; 

(3) for i £ [0,T] and u £ T>i, the map t \— > FfU is a weak entropy solution 
to CCED fort£ [0,T-t\; 



(4) if S is the SRS generated by (II. 2p . then for t £ [0,T[ and u £ V>i, 

1 



lim — 

t^o t 



F t u- (S t u + tG(u)) 



0; 



(5) fort,s £ [0, T], u,w £ T>T-t and s <t, then 



\\F t u - F t w\\ L1 < £ • \\u - u;|| L i 
\\F t u - F s u\\ L1 < £• (l + ||u|| L i) • \t 



s ; 



(1.6) 



(6) for t £ [0, T[, u £ T>i and r £ [0, T — t] the map u(t) = F T u satisfies 



(0a) for £ £ R. lim - , 



[F#u(t)) (x)-Uf u{TU) ^,x) 



dx = 0, 
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there exists a positive C such that for all a,b with — oo < a < £ < 
b < +00 



< C-(TV(«(r);]o,6[ 



lim sup — 

^0+ v 



F#u(t)) (x) 
2 



U, 



dx 



< 



where uf,s « solves (|3.2U|) and U^ u ^ ^ solves ([3.21 1) ; 



(«(t),€) 



f7) /or i G [0, T[, i/ a Lipschitz map u: [0,T — t] 1— > is swc/i i/iai u(t) G 
Z\+t ond /or 7" G [0, T — t] 



/or e e 



I, lim - 



£-i?A 



^u(r)) 



17, 



(«(r),0 



= 0, 



([^j i/iere exists a finite measure [i T such that for all a,b with — 00 < 
a < £ < b < +00 



lim sup — 

where Iff 



b-tfA 
a+i?A 



^u(r)) (x) - <( T ), 5) (^^)|| rfx < (^r (Kb[))' 



'("MO so ^ ves P-^OP an ^ (u(t) i) so ^ ves P-2ip . i/ien u(i) = 
^t«(0). 

Moreover, i//i,/2 bot/i satisfy (F) and G\,G2 both satisfy (G), then, de- 
noting by F % the process generated by fi and Gi, for all t G [0, T] and u G T>q 



F}u-F 2 t u 



L 1 



< ^.p/i-D/allco 



+ £ • ||Gi - G ! 2|| c o(^ o .i J i( K . K n)) • i- 



(1.7) 



For the definition and properties of the SRS, refer to .6]. Point 21 char- 
acterizes the tangent vector to 1 1— ► i^tt in the sense of (5J § 5]. It is through 
this characterization that the integral inequalities O and [7] are proved. The 
proof of Theorem 11.21 follows from the results presented in Section [3] below. 

The first part of next section is devoted to the application of the above 
result to the Euler system for a radiating gas and to Rosenau regularization 
of the Chapman-Enskog expansion. The framework of local sources is re- 
covered in the subsequent paragraph and, finally, we quickly comprise also 
the case of a non autonomous source. 

Remark that in the estimate f 1 1 . 6 [> the presence of the term ||u|| L i is 
mandatory, as the example dtu = u shows. Indeed, the domains Us is 
unbounded in L 1 . Moreover, note that the analogous estimate in [2"1 [8"1 [TT] 
should be understood with a time Lipschitz constant dependent on the L 1 
norm of the initial datum. 

We stress that the estimate (I6bj) is sharper than [1, formula (5.18)] thanks 
to the finite total variation of the source term, ensured by (G). 
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2 Applications and Extensions 
2.1 Euler System for a Radiating Gas 

The following model for a radiating polytropic gas was considered in |24} 
Chapter XXII, § 6], see also [TU formula (1.2)]: 

dtp + d x (p v) = 

d t (pv) +d x (pv 2 + p) = 

d t (pe + \pv 2 ^j +d x (v [pe + \pv 2 + p) + q) = 
-dl x q + aq + bd x $ A = 

Here, as usual, p is the gas density, v its speed, e the internal energy, p the 
pressure, i9 = e/c„ the temperature and q is the radiative heat flux. The 
system is closed by means of the equation of state and specifying the values 
of the characteristic constants a and b. 

Solving the latter equation in q we have q = —-^=Q a * ^^r"!? 4 ^ , where 

Qa(x) = \ exp (— y^lxQ and we are lead to consider the system 
d t p + d x (pv) = 

< dt(pv) +d x (pv 2 +p) = (2.1) 




It is well known that Euler system satisfies (F). Condition (G) holds by 
Proposition II. H Hence, Theorem 11.21 applies and we obtain the local in 
time well posedness of (I2.ip . Note that this result also ensures the local 
Lipschitz dependence of the solutions to (|2.ip from the parameters a and b. 

2.2 Rosenau Regularization of the Chapman-Enskog Expan- 
sion 

In his classical work [20J , Rosenau proposed a system of balance laws that 
provides a regularized version of the Chapman-Enskog expansion for hydro- 
dynamics in a linearized framework. The ID version is the following: 

/■ 

d t p + d x v = 
, d t v + d x p = p* * d xx v 

dt (f + d x v = X, * d 2 J 

where p is the fluid density, v is its speed and ~d is the temperature, p*, 
respectively A*, is a convolution kernel related to viscosity, respectively to 
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thermal conductivity. This linear system motivated analytical results, see 
for instance \15\ \T7\ [21] , mostly related to the quasilinear scalar equation 

dtu + d x ( \u 2 J = — u + Q *u 



since the source term — u + Q * u is equal to Q * d xx u, provided Q{x) = 
2 exp (— \x\) . Therefore, it is natural to consider the following Euler system 
with Rosenau-type sources 



i2 „. 

(2.2) 



d t p + d x (pv) =0 
d t {pv) + d x {pv 2 + p) = /i, * d\ 

d t [pe + ^pv 2 ^ + a x (v (pe + \pv 2 +p s j+q S j = A* * 

Rosenau kernels, see [201 formulas (4a) and (6)] read 

P*{ x ) = „ - exp (— bl/e) and XJx) = exp (— b|/e) 

for suitable positive parameters A, m, s, e. With the above choices, the 
sources in the last two equations in (|2.2[) can be rewritten as 

p**d 2 v = (——v + p* * v) and A**dL# = f — — i? + A* * •& 

By Proposition ll.il system (|2.2I) falls within the scope of Theorem 1 1.2 1 Thus, 
we prove the local in time well posedness of (j2.2j) as well as the local Lipschitz 
dependence of the solutions to (|2.ip from the parameters p, A, m, s, e. 

2.3 Local Inhomogeneous Source 

Theorem 11.21 can be applied also in the standard case of a local source. 
Indeed, it is immediate to see that (G) is implied by the following condi- 
tions (gl) and (g2). 

Proposition 2.1 Let j:1xOh W 1 be such that 

(gl) there exists an L\ > such that for u\,U2 € £1, \\g{x, 112) — g{x, iii) < 
L\ • ||«2 - '"ill; 

(g2) there exists a finite measure p onM such that for and x\,X2 £l 

with x\ < x 2 , \\g{x2-,u) - g{xi,u)\\ < p (]xi,x 2 ]). 

and assume that f satisfies (F). Then, setting {G{u)j (x) = g(x,u), Theo- 
rem \l.S\ applies. 

Note that the integral estimates [6] and \7\ in Theorem 11.21 ensure that the 
solution constructed here coincide with those in Similarly, the charac- 
terization S] of the tangent vector imply that the present solutions coincide 
with those in [21. 
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2.4 The Non Autonomous Case 

Theorem 11.21 can be extended to the non autonomous balance law 

d t u + d x f(u) = G(t,u) 
provided / satisfies (F), G: [0,T o ] x L 1 h L 1 satisfies 

(C) For positive 5 ,T , the map G: [0,T o ] x U So ■-> L^R.R 71 ) admits suit- 
able positive L\,Ij2,L3 such that for all u,w 6 ^ D and for all t,s £ 
[0,T o ] 

- G(s,io)|| Ll < Li • (||«- w|| L i + |i-s|) 
TV(G(t,it)) < L 2 -TV(n) + L 3 . 

Indeed, let A be an upper bound for all moduli of characteristic speeds, 
i.e. A > sup|| u || <( 5 o maxj = i r .. )n |Aj(n)|, and define 

f(u,w) = (/(«), A u/) G(u,w) = (g(J r w,u) ,x [Q1] 

here, \, , is the characteristic function of the real interval [0,11. Then, 

' A [o,i] L ' J ' 

/ satisfies (F) and G satisfies (G), so that Theorem 11.21 applies and the 
balance law dt(u,w) + d x f(u,w) = G(u,w) generates the operator F. The 
Cauchy problem 

f d t u + d x f(u) = G(t,u) 
) u(0,x) = u (x) 

is solved by t i— > F t (u , 0) where F t (u , 0) is given by the first n component 
of F t (u o ,0). 

Again, the integral estimates [6] and [7] in Theorem 11.21 ensure that the 
solutions constructed here coincide with those in |llj . 

3 Technical Proofs 

Proof of Proposition The Lipschitz property is immediate. To prove 
the bound on the total variation, call Lip (h) the Lipschitz constant of h. 
Then, it is sufficient to compute: 

TV{Q*h(u)) < supV / \\Q(y)\\ \\h(u(xi-y)) -h(u(x^i-y)) dy 
< Lip (h) sup / \\Q(y)\\ V] \\u{xi - y) - u(x i _ 1 - y)|| dy 

JR ; 
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< Lip(/i)sup / \\Q(y)\\TV(u)dy 

< Lip(fc)||Q|| L1 TV(«) 

□ 



3.1 Convective part 

Let Ai(u), \2(u), . . . ,Xn(u) be the n real distinct eigenvalues of Df(u), in- 
dexed so that Xj < Aj+i for all j and u. The j-th right eigenvector is rj(u) 
and we assume that ||r.,(0)|| = 1. 

Let a i— > Rj(a)(u) and a i— > Sj(cr)(u) be respectively the rarefaction and 
the shock curve exiting u. If the j-th field is linearly degenerate, then the 
parameter a above is the arc-length. In the genuinely nonlinear case, see [6j 
Definition 5.2], we choose a so that for a suitable constant kj > 

—Xj (iJ»(u)) = k 3 and ^(0)(0) = r,(0) 

(Sj(*)(u)) = kj, and ^(0)(0) = r,(0). 

The above choices were introduced in [2J § 2] , see also [61 [7] . 
Introduce the j-Lax curve 



Rj((r)(u) if er > 
Sj(a)(it) if < 



and define the map 

*(ct)(m~) = Vn(o-n) o . . . o -0l(cri)(u~) . 

By [51 § 5.3], given any two states u",u + € f2 sufficiently close to 0, there 
exists a vector (a±, . . . ,a n ) = E(u~,u + ) such that u + = ^>(a)(u~). 
Similarly, let S be defined by 

u + = S(a)(u~) = S n (a n ) o ...o Si(ai)(u~) 

as the gluing of the Rankine - Hugoniot curves. 

For a sufficiently small <5 Q , let u € U$ o be piecewise constant with finitely 
may jumps sited in a finite set of points denoted by I(u). Let o x ^ be the 
strength of the i-th wave in the solution of the Riemann problem for (jl.2p 
with data u{x—) and u{x+). i.e. (cr Xi x, . . . , cr x ,n) = E (u(x—),u(x+)j. Ob- 
viously if x 1(u) then = 0, for all i = 1, . . . , n. As in [6l § 7.7], A(u) 
denotes the set of approaching waves in u: 



A(u) 



{(x,i),(y,j)) e (l(u) x{l,...,n}) 2 : 
x < y and either i > j or % = j, the i-th field 
is genuinely non linear, min {(r x ,i, a y,j} < 
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while the linear and the interaction potential are 

n 

V(u)= ^ ^2\(Tx,i\, Q(u) = ^ \ a x,i^yj\- 

x£/(ti) i=l ((x,i),(y,j))£A(u) 

Moreover, let 

T(u) = V(u) + C Q(u) (3.1) 

where Co > is the constant appearing in the functional of the wave-front 
tracking algorithm, see [HJ Proposition 7.1]. Finally we define 

V$ = cl ju G L 1 (R, IR n ) : u is piecewise constant and T(u) < <5 j (3.2) 

where the closure is in the strong L 1 -topology. We remark for later use that 
there exists a positive constant c = c(5 ) with c G ]0, 1[, such that for all 
Se[0,5 o ] 

Proposition 3.1 Let f satisfy (F). Then, there exists a positive and suit- 
ably small 5 such that (jl.2p generates a Standard Riemann Semigroup 
(SRS), with Lipschitz constant L, defined on the domain T>s, for all 5 G 
}0,S o [. 

We refer to [HI Chapters 7 and 8] for the proof of the above result as well 
as for the definition and further properties of the SRS. 

Lemma 3.2 Let f satisfy (F), Q be a sufficiently small neighborhood of the 
origin; a,b £ W 1 and s G [0, +oo[ be sufficiently small. Choose u~,v~ G f2 
and define 

u + = u~ + sa , v + = v~ + sb . 

Then, there exist o~ , o~ + such that v~ = fy(a~)(u~) and v + = ^>{o~ + )(u + ). 
Moreover, 

n I n 

^2\4-^\<0{l)- i\\ a -b\\ + J2\^ )•«■ (3-3) 

i=i y i=i 

An entirely analogous result holds with replaced by S, i.e. v~ = S(o~~)(u~) 
and v + = S(a + )(u + ). 

The proof is an extension of [2j Lemma 2.1] and, hence, omitted. 



9 



3.2 Source part 

Concerning the source term we have the following results. 

Proposition 3.3 Let G satisfy (G). Then, for any 5 E ]0, <5 [ and u £ Us, 
the Cauchy Problem (jl.3p with initial data u admits a solution T,fU defined 

for t £ [0,T] where T = min j s & £~+ L3 , L ^ +1 } ■ Moreover the trajectory T^u 
has the following properties for all t 6 [0, T] and u, v G U$: 

||Et«— Stto|| L i < \\u - w\\ L i e Llt 

TV(£ 4 u) < 6 + (5 L 2 + L 3 ) t. 

The existence of S^u is a standard application of Banach Fixed Point 
Theorem and the estimates follow from Gronwall Lemma. 

Since the computations on the convective part is mainly done on piece- 
wise constant approximate solutions, we need to approximate the source 
term with piecewise constant functions. 

Let PC(M; W 1 ) be the set of piecewise constant functions in L . For any 
N e N, define the operator IIj V :L 1 (M;M n ) ^ PC(l;l re ) by 



-1+N 2 M+l)/N 



-1-N 2 

Lemma 3.4 LTjv is a linear operator with norm 1. Moreover, TV (Unu) < 
2TV(u) and for all u G L 1 (R; W 1 ) D BV(M; 1"), U N u -> u in L x (R;R n ). 

Proof. Linearity and the estimate on the norm are immediate. 

tv(il nU ) = ||(n w «)(-jv)|| + ||(n^«)(jv)|| 

N 2 -l 

+ \\(Ti N u){(k + l)/N) -(U N u)(k/N) 

k=-N 2 

r—N r-N 

< N / \\u(Z)\\dZ + N / \\u(0\\dt 

J-N-l/N JN-l/N 
N 2 -l 



k=-N 2 



r(k+l)/N rk/N 
Jk/N J(k-1)/N 

< TV(u;]-oo,-N[)+TV(u;)N-l/N,+ooty 

^ 2 -l n(k+l)/N 
+ N E / \\u(()-u(t;-l/N)\\dt; 

k= _ N 2 Jk/N 

< TV («; ]— 00, — JV[) + TV (u; 1 N — 1/N, +00 f) 
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N 2 -l 

+ Y, Tv(u;](k-1)/N,(k + 1)/N\) 

k=-N 2 

< 2TV(«). 



Concerning the pointwise convergence Hn — > Id: 

/—N—l/N p+oo 
\\u(0\\<%+ / ||u(£)| 
-oo JN 



N 2 -l 



" r(k+l)/N r(k+l)/N 

+N V / / ||u(0-«(a;)||dfda; 

/ — N—l/N r+oc i 

/ ||«(OH + ^tv(«) 
-oo JiV iV 



< 



So that limjv-»+oo ||Hjv"W — ^||x,i =0. □ 

Corollary 3.5 Lei G satisfy (G). Then, for any N, also Ily ° G satis- 
fies (G) with the same Lipschitz constant L\ and with L2,L 3 replaced by 
2L2 and 2L3. 

To simplify the operator splitting algorithm, we substitute the semi- 
group generated by (jl.3p with the Euler approximation Pt'.lts l— * L 1 (IR;M n ) 
defined by 

P t u = u + t G(u) . 

It is immediate to prove that Pt is L 1 -Lipschitz with constant 1 + tL\. We 
often use below the estimate 1+L\t < e Llt . On the other hand, observe that, 
despite the notation (useful in the sequel) , P does not satisfy the semigroup 
composition law. Indeed we can only say 

\\P s P t u - P s+i n|| Ll < 0(1) • s ■ t • (1 + \\u\\ h i) . 
3.3 Operator Splitting 

Lemma 3.6 Let G attain values in PC(M;lR n ). Then, there exist positive 
5 and T Q such that for all s € [0, T D ] and for all piecewise constant u £ T>g 

V(P s u) < V(u) + O(l) s [L 3 + V(u)] (3.4) 
Q(P s u) < Q(u) + 0(1) s\L 3 + V 2 (u)] (3.5) 
T(P s n) < T(u) + 0(l)s[L 3 + V(u)] (3.6) 
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Proof. Let u' = P s u, so that u' = u + sG(u). Call a' xi (resp. a X; i) the 
size of the i-wave in u' (resp. u) at the point x, observe that a' x i (resp. cr Xj j) 
vanishes whenever x X(w') (resp. x I{u)). Then, by Lemma f3.2l and (G) 



V(u')-V(u) < E E 

x£l(u')Ul(u) i=l 



x£l(u')Ul(u) \ 1=1 



< 0(1) • s • f V(u) + TV (G« 

< 0{l)-s-(L 3 + V(u)) . 

To derive (|3.5p . we observe that if {(x,i), (y,j)) £ A(u')\A(u), then either 
o~' x fix^i < or (y'yjOy^j < 0. Suppose, for instance, that a' xi a x ^ < (the 
other case being similar) then by (|3.3f) we obtain 



+ \<7x,i 
and therefore 



< 0(1) s ( ^l^il + || AG ( U )( 5 



E 

((*,»),(» j))e^(«')\-4(«) 



<Txj<Tyj | < O(l) s (V(u) + L 3 ) . (3.7) 



Applying (|3.7|) . (|3.4|) and since s, V(tt), y(u') << 1 we finally get 



A(u') 



y-j 



s E 

A{u') 

< 



a 



v-j 



+ Wx.; 



a yJ a vj 



' A{u') 



0(1) s { [^(«) + L 3 ] V(u') + V(u) [V(u) + L 3 ] } 

+ Yl \ a x,i(ry,j\ + E l "^' 
A(u')\A(u) A(u) 

< 0(1) ■ s ■ ( V(u) 2 + L 3 ) + Q(u) . 



i a y,j I 



Since L 3 is a possibly null constant which depends only on the system, we 
can say L 3 + V(u) = 0(1). Finally, the latter estimate follows combining 
the previous results. □ 
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Corollary 3.7 Let 5 G ]0, S [ and assume that G satisfies (G). Then we 
have 

P S F>8 Q ^S+0{l)s(2L 3 +S) and PsUs Q Us +s (l 2 6+L 3 ) ■ 

In particular take a constant C > 0(1) (2L3 + <5 Q ) , a number 5 G ] 0, <5 D [ and 
time T G ]0, T a ] suc/i i/iaf 5 + C T < 5 a then for any t G [0, T\, s G [0,T — t] 

Ps^S+tC Q ^«5+C(t+s) (3-8) 

Proof. Fix it G and an approximating sequence of piecewise constant 
function with T(ufe) < <5. The previous Lemma shows that 

T(u A + s(IItfoG)(ufc)) < T(u k ) + 0(l)s[2L 3 + V(u k )] 

< 5 + 0(l)s[2L 3 + 8] . 

But Mfc + sLTiv G(u k ) converges to P s u as k,N — > +00 and so P s u G 
^ > (5+0(i)s(2L 3 +(5)- The proofs of the other inclusions are straightforward. □ 

Corollary 13.71 allows to define the domains appearing in Theorem 11.21 as 

V t = V 5+ ct VtG[0,T]. 

Let h G N and define 

F?u = St-UP>°Ss) h u t€ + ( 3 - 9 ) 

In other words, in any interval ]hs, (h + we apply the semigroup 5. In 
turn, at the times t = hs, P s is applied. 

If a time T and a 5 G ]0, <5 [ are chosen as in Corollary 13.71 then Ffu is 
defined up to the time T — t for any u G T>i and i G [0, T\. 

Observe that for t' , t" G [0, T] with t' + t" G [0, T], the following inclusion 
holds: 

F$V t » C . (3.10) 

Note that P t s is L 1 -Lipschitz with Lipschitz constant bounded by L-L t ^ s e Llt , 
with L as in Proposition 13.11 and L\ as in (G). The following theorem 
shows that the Lipschitz constant of Ff actually is bounded from above by 
a quantity independent from s. 

Theorem 3.8 Let f satisfy (F) and G satisfy (G). If 5 Q is chosen suffi- 
ciently small and 5, T are chosen as in Corollary \ 'J. 7\ then there exists a 
constant C such that for all i G [0, T], u, w G T>i and t G [0,T — t\, we have 

\\F t s u - F t s w\\ L1 <C-\\u- w\\ L1 . (3.11) 
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Proof. The proof can be carried out following essentially the same line used 
in [H Theorem 4.2]. Therefore we only outline it. The key ingredient is 
the Liu & Yang functional which, unfortunately, is defined only on piece- 
wise constant ^-approximations of the trajectories of the convective part 
(see [6]). Therefore we first suppose that the source G(u) attains values in 
PC (R,M n ). For this kind of sources one can prove (|3.1ip following almost 
exactly the proof of [21 Theorem 4.2]. The only difference is the global nature 
of our source term. But the non locality of G{u) can easily be tackled using 
Lemma 13.21 integrating and summing up the pointwise estimates obtained 
in [21 Theorem 4.2] and finally using (G) similarly to what we have done in 
a detailed way in the proof of Lemma 13.61 

Once we obtained (13. lip for piecewise constant source terms, we apply 
it to the source term IIjv o G{u). If we denote by F^ ,N u the trajectory 
defined by (|3.9I) with II at ° G{u) in place of G(u), it is quite easy to see 
that Lemma [3.41 implies the strong convergence of F^ ,N u to Ffu proving the 
theorem for a general source G. □ 



We need below the following estimates concerning the dependence of F s 
on time. 



Lemma 3.9 Let T and 6 be as in corollary \3. 7| Then, for all u &T>i and 
allte [0,T-t], 



(PsTu 



\ F t u - u \\-Li 



{Psfu - U 



< 0(i)-(i + HIlO 

< 0(i)-(i + NIli) 

< 0(l)-t-(l + ||n|| L1 ) 

< 0(1) ■ k ■ s ■ (l + |M| L i) 



(3.12) 
(3.13) 
(3.14) 
(3.15) 



Proof. Consider first (|3.13j) with k = 1. By (G) and the properties of S 

\\P s u\\ h i = ||u + sG(it)-sG(0) + sG(0)|| L1 

< ||u|| L i + sLi||u|| L i + s ||G(0)|| Ll 

< e sLl (cs + ||u|| L i) 

H-SWIIl 1 < \\S s u — «|| L i + IMIli 

< L s + |M| L i 

< e Ls (cs + ||u|| L i) 

where c = max{L, ||G(0)|| Ll } and L is as in Proposition 13.11 Proceed now 
by induction on k and, for fcs£ [0, T — t\, 

||F fc >|| L1 < e^ L+L ^(2cks + \\u\\ Ijl ) 



< e ^ L+L ^ ks (2cT+ ||«|| Ll ) 

< O(1)-(1 + Hl0 • 
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Therefore, for i £ [0, T — t [ and k = [|] , we have 

< Ls + 0(1)- (1+Hn) 

< o(i)-(i + Ni L i) . 

The estimate ()3.13p is obtained similarly. Passing to the Lipschitz estimates, 
for (k + l)s G [0, T — t[, by (l3TT2|) 



L 1 



|^>-n>|| L i + ||^ 



u — u 



L 1 



< ||5 s F fe > + sG(S s F^u) - F fc >|| Ll + 



M — u 



L 1 



< S • (L + Lt \\S s F^u\\ L1 + ||G(0)|| L1 J + \\F% s u - u|| Ll 

< 0(l)-a-(l + |H| Ll ) + ||F^u-u|| Ll . 



By induction, ||i^tt — ^|| L i < C(l) • /c • s • (l + ||u|| L i). 
We are left with the case - g" N. If t G [0, s[ we have 



|F/n - u|| Ll = ||5 t « - «|| L i < Oil) ■ t , 



while if t > s, so that k = [-] > 1, 

- it|| L i < Ls + 



Ff u-u 



L 1 



< Ls + 0(1) -k-s- (l + ||m|| l i) 

< 0(l)-t-(l + N| L i) • 



The proof of (|3.15|) is entirely similar. 



□ 



The next step consists in showing the convergence of F s as s tends 
to zero. This result will be obtained with the help of the commutation 
relation (jl.5p which we will show to be true for St and Pp We need the 
following result that is an easy consequence of {2j Remark 4.1]. 

Proposition 3.10 Take u,v,u G T>t, then, for any t\ < t<i one has the 

estimate 



\S t , 2 w - S t2 u - uj\\ l1 < L ■ \\S tl w — S tl u 



L 1 



(3.16) 



+ 0(l)-(t2-ti)-TV(w). 
Proposition 13. 101 implies the following commutation result: 



Theorem 3.11 Let 5 and T be in Corollary \371[ For any t G [0,T[, u G T>i 
and t G [0, T — i] , we have the estimate 



\\StPtu-P t S t u\\v<0(l)-t 2 



(3.17) 



15 



Proof. Since u + tG(u), tG(Stu) and u all belong to T>t, we can apply (|3.16|) 
with uu = u + t G(u) and u = tG (Stu) to obtain: 



S t P t u - P t S t u\\v = S t [u + t G(u)} - Stu-tG (S t u) 

< L ■ S [u + tG(u)} - S u-tG(S t u) 

+ 0(1) -i -TV [tG(S t u)] 

< L ■ t ■ \\G(u) - G (S t u)\\ jA + Oil) ■ t 2 

i • 



< L 2 ■ Li ■ t 2 + Oil) ■ t 2 



□ 



Now we show that (|3.17|) and the uniform Lipschitz property (|3,lip of 
the approximations imply the existence of a "tangent vector" and the strong 
convergence of the approximations. We will show that these two conditions 
are enough and that there is no need to use again the almost decreasing 
functional as was done in Lemma 5.1]. 

Proposition 3.12 Let 5 and T be as in Corollary \'J. 1\ For any t G [0, T[, 

u G T>i, t G [0, T — t\ and s, s' € ] 0, i 2 ] , we have 



\F t s u - StPtuL, < O(l). (l + H L i) -t 2 



\F t s u-P t S t u 
F t s u-F?u 



L 1 



< Oil)- fl+lklU -t 2 



(3.18) 



± < 0(l)-(l + ||n|| Ll )-t 2 



Proof. We prove only the first inequality in (|3.18[) , the other two inequalities 
being consequences of this one and of Theorem 13.111 For any integer k £ 
[1, T/s] define 



Pk{s) 



sup 

u£T> T _ ks 



(P s ) k S s u-S s (P s ) k u 



Now, for u G 2V-(fc+i)s (and hence P s u G T>T-ks), we can compute 



< 



'P s ) k+1 S s u-S s (P s ) k+1 u 
[P s ) k P s S s u-iP s ) k S s P s u 



< 



+ 



'PsYSsPsU-SJPsTPsU 



< e 1Ll \\P s S s u- S s P s u\\ Ll + sup (P s yS s w - S s (P s ) k w 

< e TLl Pl (s) + p k (s). 

And hence Pk+i{s) < & TLl p\(s) + Pk(s) that, by induction, gives 

Pk(s) < e TLl k Pl (s). 



L 1 



16 



Now define 



Pk(s) 



sup 

u£V T _ ka 



F^ s u-S ks {Ps) k u 



L 1 



Again we can compute for u £ T> 



T-(k+l)i 



%+i) s u - S(k + x) a (P s ) 



fe+1 



u 



< 



L 1 



< \\F£ s P s S s u-F£ s S s P s u\\ L1 + 



Pks^sPsU — Sks(P s ) k S s P. 



u 



L 1 



+ 



Sks(Ps) S S P S U SksS s (P s ) Ps 



L 1 



< Cpx{a) + p k (s) + Lj(P s ) k S s P s u - S s (P s ) k P 

< Cpi (s) + p k (s) + Lp k (s) 

< (c + Le L ' T k} Pl (s) + p k (s). 



u 



L 1 



Therefore we have Pk+i(s) < ^£ + Le LlT fcJ + Pfc(s) which gives, to- 
gether to pi(s) = pi(s), by induction 

Phis) < (£ + Le LlT k) k Pl (s). 



Fix now t £ [0, T - t ] , take se]0,i 2 [ and define k= [|] . We have for all 



T-ks 



F t s u - S t (P s ) k u < 0(l)-(t-ks) + 



L 1 



A t 



L 1 



< 0(1) • s + 0(1) • (ts + t 2 

< 0(l)-t 2 , 

where the last inequality is a consequence to the fact that Pl ^ is bounded 
because of (|3,17p . 

We are left to prove that 



(P s ) k u - P t u | < 0(1) ■ (1 + ||u|| Ll ) • t 2 . 



By HEED, 



(P s ) fc+1 u-P ( 



t/ 



L 1 



P s (P s ) fe u-P fes u- S G(u) 



L 1 



(P s ) fc u - P fcs « + s G ((P s ) fc n) - s G(u) 



L 1 



< 



(P s ) fc u - P fcs « + 0(1) • (1 + |M| L1 ) • k s 2 
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We thus recursively obtain 



[PsV 

Therefore 

{P s fu - P t U 



PksU 



< 



x < C(l) • (1 + ||^|| L1 ) • A; 2 s 2 . 



{P s ) k u - P % u 



+ 



P %s u - P t u 



< 0(1)- (l + |H| L i) ■k 2 s 2 + 0(l)-s 



< 0(1)-(1 + IM 



Li 



■t 2 . 



□ 



Now we prove the convergence of the approximations and the character- 
ization of the tangent vector, i.e. [H in Theorem 11.21 

Theorem 3.13 Let 5 and T be as in Corollary \3.1\ For any i G [0,T[, 
« £ Dj and t G [0, T — 1\ the sequence Ffu converges in L 1 as s — ► to a 
limit trajectory F^u which satisfies the tangency conditions 



\\F t u-S t P t u\\v < 0(1)- (l + |H| L i) -t 2 
\\F t u-P t S t u\\^ < O(l). (1 + |H| L1 ) -t 2 



(3.19) 



Proof. Because of (I3.18p . we need only to show that s — > Ffu is a Cauchy 
sequence in L 1 as s — > 0. Fix e > arbitrary. Then choose = to < t\ < 
. . . < ijv— i < tiv = t so that tj — ti-i < e for i = 1, . . . , N. Then observe 
that Definition (|3.9p . Theorem 13.81 and f)3. 12|) imply that F s satisfies an 
approximated semigroup condition: 



F^u-F^u 



L1 = o(i).(i + H L1 ). s . 

Therefore, for any < s, s' < miiij = i..jv \{U — ij_i) 2 }, we can compute 

N 



Ffu-Ffu 



L 1 



i=i 



L 1 



< F s F s 'u — F s F s F s ' 



u 



i=l 



L 1 



+ o(i) (i + |M| L1 )iVs 



iV 



< £EK U "^^ U Ll + o(i)(i + N| L1 )iv s 

i=l 

/ ^ L z L z — 1 H— 1 L a H— 1 H— 1 



i=l 
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+ 0(1) (l + H L i) N(s + s') 



N 



_ 0(iWi + |M 



8=1 



< 0(1) (l + H| L i) {et + N(s + s')) . 
And, finally, as s, s' — ► we get 



lim sup 

s,s' -s-0 



F t s u-F t s u 



x <0(l) (l + N| L i) et 



which proves the Theorem because of the arbitrariness of e. 



□ 



The limit trajectory thus obtained satisfies © in Theorem II, 2 \ as can be 
seen passing to the limit s — > in (13.101) and in the approximate semigroup 
condition 



Ff^u - F? i+t2 u 



i < O(l) - (1 + ||n|| Ll ) - S . 



Taking the same limit in (|3.1ip . we prove the former inequality in (II. 6p . To 
prove the latter estimate, observe that by (|3.14p 

||i*iv-«|| L i < 0(1) • (1 + HlO 
while, for t<i > t\, the semigroup property implies 

\\F t2 u - -F^Hli < \\F t2 - tl F tl u - F tl u\\ Ll 

< 0{1)- (l+||F tl «|| Ll ) -{t2-h) 

< 0(l)-(l + H Ll )-(i 2 -ii). 

Assertion follows from (|3.19p . 

We pass now to ([3]). The trajectory t \— > F t u is a weak entropic solu- 
tion of (jl.ip . This can be proved using the properties of the approximate 
solutions constructed above, as in (HEJ[Tl]. Here, we prefer to exploit the 
tangent vector provided by theorem 13. 131 



Corollary 3.14 Let 5 and T be as in Corollary \3.1\ For any t £ [0,T[, 
u £ T>i and t G [0, T — t ] the trajectory 1 1— > Ftu is a weak entropic solution 



For the definition of weak entropic solutions of a balance law, refer to [12j . 
[E2 (2.16) and (2.19)] or [11, § 6]. 

Proof of Corollary \3.14\ We show below only the entropy inequality, since 
the proof that 1 1— > Ftu is a weak solution is entirely similar. 



19 



Observe that, by (13131) 

F t u = S t u + 0{1) ■ (1 + ||u|| L i) -t 
F t u = S t u + tG(u) + <D(l)- (l + H L i) -t 



2 



in L . Let (77, q) be an entropy-entropy flux pair and tp € be a non 
negative test function. Fix a positive e and denote Ii = [ie, (i + l)e[ x M. 
for i G N. By the properties of 5, dtrj(Stu) + d x q(Stu) < in the sense of 
distribution, using the Divergence Theorem we get 

n(r](F t u)dt<p + q(F t u)d x (p) dxdt = 
= YlJJ [v( F tu)d t (p + q(F t u)d x (p) dx dt 

i * 

= ^ y J J (r)(S t -ieF ie u)dt(p + q(St- i£ F i£ u)d x ip) dx dtj 

+ 0(l)(l + |M| L1 )e 
> Yl (J w (jl{S E F iE u) (p ((i + l)e, x) - r)(F i£ u) <p(ie, xfj dx\ 

+ 0(l)(l + H Ll )e 

= Yl (j R {v(. S eFieU) ~ v( F (i+l)eU)) <f ((* + l)e, x) dx 

+ <5(l)(l + IHI L i)e 

= E (X ^ (5£jFfeU) " ^FeFieU)) V ((i + 1)E, x) dx 

+ 0(l)(l + H Ll )e 
= E / (^^-^(^FfeU + EGCi^u))) y>((i + l)e,s) 

+ 0(l)(l + |H| Ll )e 
= J2 f (-Dv(Fieu)eG(F ie u)<p((i + l)e,x)^ dx 

+ <9(l)(l + H Ll )e 
= ^J£[-DTi(F t u)G(F t u)<p(t,x)) dxdt 

+ 0(l)(l + |H| L1 )e 
= / / -Drf(F t u)G(F t u)<p(t,x)dxdt + 0(1) (l + ||u|| L i) e 
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By the arbitrariness of e, we conclude with the distributional inequality 
d t r]{F t u) + d x q(F t u) - D V (F t u) G(F t u) < . 

□ 

Now we show a result on the dependence of the solution with respect to 
the source term. 

Theorem 3.15 Let f lf f 2 satisfy (F) andG x ,G 2 satisfy (G). Call F 1 ,F 2 
the corresponding semigroups and assume they are defined on a common 
family of domains T>t- Then, for any t G [0, T[, u G Vf and t G [0, T — t], 
the Lipschitz estimate (|1.7j) holds. 

Moreover, fix a flux f satisfying (F) and a sequence of source terms 
Gk satisfying (G). If G^ converges pointwise to G, then the corresponding 
semigroups F k converge pointwise to the semigroup F generated by G. 



Proof. We apply the well known integral estimate, see [H Theorem 2.9]: 



Flu-F?u 



< £ I liminf 



F$ +T u-F*F?u 



o 



dr 



By (|3.19p and using the stability result [H Corollary 2.5], with obvious 
notation we have 



lim inf — 

< lim inf — 

< lim inf — 



FlF 2 T u-FlF 2 T u 



< 



L 1 



S$Ffu + 0G 2 (F>) - S^F'u - 0Gi(*>) 



+ G 2 {F$u) - Gx(F 2 u 

< 0(1) ■ (J|-D/l - £>/2|| C 0(n,R™ x ") + 11^1 ~~ ^2 1 1 QO Qi &o jL 1 (M;M n )) 

Concerning the pointwise convergence, note that 

F t k u -F t u <cf \\G(F T u) - G k (F T u)\\ x dr 
Ll Jo 



and the proof is concluded trough Lebesgue convergence theorem. 



□ 
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3.4 Integral Characterization 

Following § 9.2], P § 5.2] and [3, Definition 15.1], let v G V T . For all 
define uf v ^ as the solution to the homogeneous Riemann problem 



d t w + d x f(w) = 
w(0, x) = 



lim v(y) if x < £ 
lim u(y) if x > £ . 



(3.20) 



Define E/t^ as the broad solution (see [U § 3.1]) to the Cauchy problem 



(3.21) 



j + Df (v(0) d x w = G{v) 
I w(0, x) = v(x) , 

so that denoting if = k (v(£)), rf = rj and Af = A, (u(£)), 

L^ftx) = [/ ( b 4(t,x) + ^ b ; 2 5) (t,x) 

i=i ^ ' 
n f f / \ 
O*'*) = E/ (if-G(«)(x-Af*))ff«i>. 

We are now ready to prove the first part of the characterization stated in 
Theorem 11.21 

Theorem 3.16 Let F be the map constructed in Theorem \3.13[ let A be an 
upper bound for all characteristic speeds. Then, for all t E [0, T[, u € T>i 
and all r 6 [0, T — t\ , F T u satisfies^ in Theorem 



dx < 



Proof. To obtain (|6ajl . compute: 



f+#A 



{F#u(t))(x)-U1 (#,x) 



I /-f+^A 
< - 



£-tfA 
£+tfA 



(F*u(r)) (a?) - (S*u(t)) (a?) -#(g (u(t))) (x) 



</./• 



+ 



1 

1 



£-#A 
£+tfA 

£-i?A 



(S*u(t)) (x)-t/}, i0 (i?,x; 



19 (g(u(t)))(x) 



< - 



^n(r)-^n(r)-^G( U (r)) 



L 1 
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+ 



+ 



dx 



G{u(t)))(x) 



dx 



As — > 0, the first summand above vanishes by (jH) in Theorem 11.21 the 
second by (9.16), § 9.2] and the latter one by G (w(r)) E L 1 . 
Similarly, to obtain (16b j) . we exploit [6j § 9.2]: 



< 



1 
1 

7) 
+- 



l 

< - 

- 



a+i?A 
a+i?A 



a+i?A 



+ 



+ 



{FMr)) (x)-U^ o (0,x 

[F#u(t)) (x) - [SMr)) (x) -&(G (n(r))) (x) 

{SMr)) (x) - U\ v ^ : x) + $(G (u(t))) (x) 
F, & u{t) - Svu(t) -0G («(t)) 

(5^(r)) (a?) - [/^(^x 

$(g {u(t)))(x)-U^,x) 



dx 
dx 



L 1 



a+tfA 
fe-i?A 



dx 



a+i9A 



dx 



1 

< - 

- 



L 1 



F^u(r) - 5,?u(r) - G (u(r)) 
+ 0(1) (TV («(r);]o,6 

*(g (u(t)))(z)- 17^(0,*) 



1 



a+i?A 



dx 



As i? — ► 0, the first summand above vanishes by (J2J) in Theorem 11.21 while 
the latter vanishes as — > 0. Indeed, using [61 Lemma 2.3], 




1 

< - 

- 



a+i?A 



0[g{u(t)))(x)-U^ o (0,x) 



dx < 



EfMGM) (x))r|^d S 



>,0 



i=l 

n 



■d 



^ / ^ ) -(G(n(r)))(x-A^ ) S )r^d S 
i=i J ° v y 

G(u(t)))(x)-(g (u(t)))(x-X^s) 



n r ti r b-®\ 




1 JO ^o+« 



dx ds 
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< 



n 

£ / Tv(G(u(r));]o,6[ 

8=1 ^° 



s (is 



< 0(l)-^.Tv(t?(«(r));]a,6[) 
completing the proof of (|6bj) . 



□ 



Theorem 3.17 Lei F be the map constructed in Theorem \3.13\ and X be an 
upper bound for all characteristic speeds. Assume i G [0, T[, u: [0, T — t\ i— ► 
2?r is Lipschitz, satisfies u(t) G and 5oi/i (I7ap and (|7bj) in Theorem \1.2\ 
hold. Then, u(t) = F t u for all t G [0, T - t]. 

We omit this proof, since it is a slight modification of [6, Part 2 of 
Theorem 9.2]. 



3.5 Consequences of the Hyperbolic Rescaling 

Given a function v. M i— ► M n and A > 0, we denote by v\ the function 
obtained by applying a dilatation to v, i.e. va(^) = v(Xx). Obviously 
v G T>$ implies v\ £ T>$. We have the following Proposition (see also 
Corollary 1]). 

Proposition 3.18 Let S: [0, T] x T> g T>$ be the semigroup generated by a 
system of conservation laws and let d: T>s x T>s i— > K + fre a distance satisfying 

d(u\, V\) = — d(u, v) for all u, v G Pj, and A > 0, 
A 

and the Gronwall estimate for a positive C: 

d(SfU, Stv) < e ct d(u, v) for all u,v G T>s and t G [0, T]. 
Then C = 0, i.e. d is non expansive with respect to S. 

Proof. If u(t, x) = (Stu) (x) is a semigroup trajectory, then also u(Xt, Ax) = 
(StU\) (x) is a semigroup trajectory. Therefore we have the equality 

(Stu) x = S±u\ for all u G T>$, t G [0, T] and A > . 
Hence we can compute for all u, v G T>$ and A > 

d(S t u,S t v) = Xd((S t u) x ,(S t v) x ) = Xd^St_u\,St_v\^ 
< Xe c *d(u\,v\) = e c *d(u,v) 
Now, letting A tend to infinity, we get the non expansive property 
d (Stu, Stv) < d (u, v) for all u,v G T>s and t G [0, T\. 

□ 
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